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The Fermi Gamma-Ray Space Telescope was launched in June 2008 and the onboard Large Area Telescope 
(LAT) has been collecting data since August of that same year. The LAT is currently being used to study 
a wide range of science topics in high-energy astrophysics, one of which is the study of high-energy cosmic 
rays. The LAT has recently demonstrated its ability to measure cosmic-ray electrons, and the Fermi LAT 
Collaboration has published a measurement of the high-energy cosmic-ray electron spectrum in the 20 GeV to 1 
TeV energy range. Some methods for performing a similar analysis to measure the cosmic-ray proton spectrum 
using the LAT will be presented with emphasis on unfolding the reconstructed proton energy. 



1. Introduction 

The Fermi Gamma-ray Space Telescope was de- 
signed for the study of many interesting topics in as- 
trophysics ranging from pulsars, active galactic nu- 
clei, gamma-ray bursts, and the indirect detection of 
dark matter. Another intriguing subject that has at- 
tracted strong interest is the use of the Fermi LAT in 
the study of cosmic rays. Early in the design stages of 
the LAT, the potential for making a cosmic-ray elec- 
tron measurement was acknowledged Recently 
the Fermi LAT Collaboration demonstrated this abil- 
ity and published a high-statistics measurement of the 
cosmic-ray electron spectrum in the energy range from 
20 GeV to 1 TeV [3[, containing about 4.5 million 
events collected over a six month period from August 
2008 to January 2009. The electron spectrum mea- 
surement can be used in the study of cosmic-ray prop- 
agation and in the constraint of the diffuse gamma-ray 
emission Q. Furthermore, it may also be possible to 
use the LAT to obtain a measurement of the cosmic- 
ray proton spectrum, and studies are being conducted 
within the LAT collaboration to explore this possibil- 
ity. We present some preliminary findings illustrating 
the prospects for such an analysis. 

1.1. The Fermi Gamma-ray Space 
Telescope 

The Fermi Gamma-ray Space Telescope was 
launched on June 11, 2008 from Kennedy Space Cen- 
ter. Two instruments are onboard the satellite. The 
main instrument is the Large Area Telescope (LAT), 
which is a pair conversion telescope sensitive to the en- 
ergy range from 20 MeV up to greater than 300 GeV. 
The second instrument onboard the satellite is the 
Gamma-ray Burst Monitor (GBM) and is designed 
to detect GRBs in the 8 keV to 40 MeV energy range. 
It consists of 12 Nal and 2 BGO detectors positioned 
around the spacecraft, and has a large field-of-view al- 
lowing it to constantly monitor the entire unoccultcd 



sky. The analyses described here pertain to the LAT 
instrument, and a more detailed description of the 
LAT detector can be found in 



2. Measuring the Cosmic-ray Proton 
Spectrum 

We are currently exploring methods to produce a 
cosmic-ray proton spectrum measurement using the 
LAT instrument. However this analysis is inherently 
more challenging than the electron spectrum measure- 
ment. The main cause for this is that the thickness 
of the calorimeter is optimized for electromagnetic 
showers rather than hadronic showers. Whereas the 
calorimeter is 8.6 Xo on-axis for EM showers, it is only 
about 0.43 interaction lengths for hadronic showers. 
As a result, most protons are only minimum ionizing 
particles (MIPs) in the LAT, and thus are not useful 
for the spectrum measurement. For the protons that 
do interact and produce a shower, the energy resolu- 
tion is much worse than for either gamma rays or elec- 
trons. In fact we have found that the reconstructed 
energy is approximately a lower limit on the true in- 
coming energy. In addition, the event reconstruction 
is more complicated than for electrons or gamma rays 
(Figure [T] and Figure [2]) . 

Here we present our characterization of the instru- 
ment response using a preliminary proton event selec- 
tion applied to events collected via the onboard diag- 
nostic filter (here referred to as the DGN filter) . This 
filter allows a prescaled sample of all hardware trigger 
types to be downlinked. An alternate and indepen- 
dent pathway that is also being utilized for event col- 
lection is the high-energy pass feature of the onboard 
GAMMA filter (to be referred to here as the High Pass 
filter) . This feature of the onboard software filter al- 
lows all events that deposit greater than 20 GeV in the 
calorimeter to be downlinked. The High Pass allows 
for higher statistics than the prescaled sample from 
the DGN, however with the DGN, there is no lower 
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Figure 1: A hadron candidate event: a large energy 
deposit per ACD tile, a small number of extra clusters 
around the main track; a large number of extra clusters 
away from the track, and a large and asymmetric shower 
profile in the calorimeter. 




Figure 2: An electron candidate event: few ACD tiles hit 
in conjunction with the track, a clean main track with 
extra clusters very close to the track (note backsplash 
from the calorimeter), and a well-defined symmetric 
shower in the calorimeter (though not fully contained). 
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Figure 3: The preliminary geometry factor for DGN filter 
collected events is plotted versus MC energy. 



from Monte Carlo simulations, with the full instru- 
ment geometry, of protons interacting within the LAT. 
A preliminary geometry factor as a function of Monte 
Carlo (MC) energy, or the true incoming particle en- 
ergy, is shown in Figure [H It peaks at approxi- 
mately 0.005 m 2 sr near 700 GeV, and falls off appre- 
ciably below 10 GeV. For comparison, the estimate 
for events collected via the High Pass using these se- 
lections peaks near 0.8 m 2 sr. The geometry factor 
for electrons f|[, which also uses the High Pass fil- 
ter, peaks around 2.8 m 2 sr and further illustrates the 
fact that most protons are minimum ionizing in the 
LAT. However it is worth noting that the preliminary 
proton geometry factor, using the High Pass, is still 
comparable to other cosmic-ray experiments, such as 
the average value of 0.15 m 2 sr reported by AMS in 
2002 0. 



2.2. Preliminary Background 
Contamination 



energy threshold, and the energy spectrum can be ex- 
tended to lower energies. These two means of data 
collection are thus complementary and will also allow 
a cross-check on the measured energy spectrum. In 
the following sections, we further discuss methods for 
measuring the cosmic-ray proton spectrum using the 
DGN pathway. For a earlier discussion using the High 
Pass feature to measure the proton energy spectrum, 
refer to [||. 

2.1. Preliminary Geometry Factor 

A preliminary proton event selection has been de- 
veloped, and the geometry factor has been calculated 



For events collected via the DGN filter, we have es- 
timated the background contamination given the pre- 
liminary selections using Fermi Monte Carlo simula- 
tions of the cosmic-ray environment the LAT encoun- 
ters, which has been modelled from previous cosmic- 
ray measurements and again includes the full geome- 
try of the detector. The dominant background compo- 
nents are found to be electrons/positrons, alphas, and 
particles with Z>2. The fraction that each particle 
type contributes to the selected event sample is plot- 
ted versus the reconstructed energy in Figure 3] Over 
most of the energy range, each background species 
contributes only a few percent or less to the selected 
event sample. 
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Figure 4: The fraction of background contamination for 
the different background components is plotted versus 
reconstructed energy using the preliminary proton 
selections applied to simulated DGN filter collected 
events. The fraction of alpha particles is plotted in black, 
electrons/positrons are plotted in red, and particles with 
Z>2 in blue. Each background component contributes 
approximately a few percent or less over most of the 
energy range. 



Figure 5: The preliminary energy response, for events 
collected via the DGN filter, is plotted as MC energy 
versus reconstructed energy. For a given bin in MC 
energy, there is a wide range of reconstructed values 
obtained. In addition, the reconstructed value can be 
seen to be approximately a lower limit on the MC energy. 
Both illustrate the challenge in measuring the energy of 
proton events. 



2.3. Energy Unfolding 

We have also studied the application of an energy 
unfolding algorithm to the reconstructed energy dis- 
tribution and the energy response of the instrument 
given the preliminary selections. The process of en- 
ergy unfolding is to calculate a distribution of the true 
incoming energies of an event sample, given the dis- 
tribution of reconstructed energies and a detector re- 
sponse matrix. It should be noted that this procedure 
does not attempt to correct energies on an event-by- 
event basis, but rather its goal is to obtain an esti- 
mated distribution for the incoming energies of the 
selected event population. A similar energy unfold- 
ing process was also performed for the LAT electron 
spectrum measurement, however given that the en- 
ergy resolution is much worse for protons, the effect of 
this process on the reconstructed proton spectrum will 
be much greater. The energy response used in the un- 
folding algorithm is plotted in FigureEJ plotted as MC 
energy versus reconstructed energy. This response has 
been calculated by applying the preliminary proton 
selections to Fermi LAT Monte Carlo proton simu- 
lations using a hard spectrum of E _1 . This enables 
more events to be simulated in the higher-energy bins 
than would a more typical cosmic-ray (CR) spectrum. 
From this response histogram, it can be seen that for 
a given MC energy bin, there is found to be a wide 
range of reconstructed energy values. In addition, it 
can also be seen that the reconstructed energy is ap- 
proximately a lower limit on the MC energy. 

To perform the unfolding procedure, we have used 
RooUnfold, a ROOT-based framework for unfolding 
Q . The unfolding algorithm takes as input the recon- 



structed energy distribution and the energy response, 
and the calculated unfolded distribution is returned. 
The unfolding procedure is applied to a reconstructed 
energy distribution taken from a MC sample indepen- 
dent from the one used to create the energy response. 
As a test of this procedure, the unfolded energy dis- 
tribution can be compared to the MC energy distri- 
bution. For this test we have used the reconstructed 
energy distribution, after preliminary selections, from 
Fermi LAT MC simulations of the cosmic-ray envi- 
ronment. It should be noted here that the unfolding 
procedure is applied here to only the proton compo- 
nent of the event sample surviving the preliminary 
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Figure 6: Comparison of the unfolded energy distribution 
(black points) to the MC energy distribution (red line) 
using preliminary proton selections applied to simulated 
DGN filter events from Fermi LAT MC cosmic-ray 
simulations. The input reconstructed energy distribution 
(blue line) is also shown. 
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Figure 7: The unfolded MC cosmic-ray proton spectrum 
(green circles) for events from the DGN filter, calculated 
by dividing the unfolded energy distribution by the 
simulated livetime and the geometry factor. For 
comparison, the proton energy spectrum reconstructed 
using the MC energy distribution is also shown (red 
triangles). 



2.5. Future Steps in Analysis 



Improved selections are being studied with the goals 
of spanning a larger energy range, achieving a higher, 
more constant efficiency, and a more constant back- 
ground contamination level as a function of energy. 
There is also a continuing effort aimed at improving 
the application of the unfolding procedure to produce 
a better agreement between the unfolded energy dis- 
tribution and the MC energy distribution. The sys- 
tematic uncertainties associated with the energy un- 
folding technique are also being estimated. Further- 
more, the background rates need to be estimated and 
subtracted from the proton candidate event rate, and 
in addition, the systematic errors from MC/Data dis- 
crepancies in the selection variables. These methods 
can then be applied to Fermi LAT data and the re- 
sulting distribution used to reconstruct a measured 
cosmic-ray proton spectrum. 



selections. The resulting unfolded energy distribution 
is shown in Figure [5] compared with the MC energy 
distribution. The unfolded distribution obtained is a 
reasonable reproduction of the MC energy distribu- 
tion with differences of approximately 20% or less. 

2.4. Reconstructed Monte Carlo CR 
Proton Spectrum Using the DGN 

The unfolded distribution shown in Figure [6] can 
be used to calculate the reconstructed MC cosmic- 
ray proton spectrum. The spectrum is calculated by 
dividing the unfolded distribution by the livetime of 
the MC simulation used and the geometry factor pre- 
viously shown. The resulting spectrum is plotted in 
Figure [7] (green circles). The reconstructed proton 
energy spectrum, calculated using the MC energy dis- 
tribution, is plotted for comparison (red triangles). 
Over most of the energy range, the reconstructed MC 
energy spectrum is reasonably well reproduced by the 
unfolded energy spectrum. 
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